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Abstract 

A  systematic  development  of  the  macroscopic  field  equations  (conser¬ 
vation  of  mass,  linear  and  angular  momentum,  energy,  and  Maxwell’s 
equations)  for  a  multiphase,  multicomponent  medium  is  presented.  It  is 
assumed  that  speeds  involved  are  much  slower  than  the  speed  of  light  and 
that  the  magnitude  of  the  electric  field  significantly  dominates  over  the 
magnetic  field  so  that  the  electroquasistatic  form  of  Maxwell’s  equations 
applies.  A  mixture  formulation  is  presented  for  each  phase  and  then  av¬ 
eraged  to  obtain  the  macroscopic  formulation.  Species  electric  fields  are 
considered,  however  it  is  assumed  that  it  is  the  total  electric  field  which 
contributes  to  the  electrically  induced  forces  and  energy.  The  relation¬ 
ships  between  species  and  bulk  phase  variables  and  the  macroscopic  and 
microscopic  variables  are  given  explicitly.  The  resulting  field  equations 
are  of  relevance  to  many  practical  applications  including,  but  not  limited 
to,  swelling  clays  (smectites),  biopolymers,  biological  membranes,  pulsed 
electrophoresis,  and  chromatography. 

Key  words:  porous  media,  mixture  theory,  electrodynamics,  averaging,  swelling 

1  Introduction 

We  attempt  to  address  the  following  issue:  in  a  composite  or  porous  medi¬ 
um,  how  do  the  electro-thermodynamic  variables  of  each  individual  constituent 
contribute  to  the  electro-thermodynamic  variables  of  the  mixture  as  a  whole? 
In  other  words,  given  information  about  individual  materials,  determine  the 
relationship  between  the  electro-thermodynamic  properties  of  individual  com¬ 
ponents  and  the  electro-thermodynamic  properties  of  the  averaged  multi-phase 

*  University  of  Colorado  at  Denver,  Center  for  Computational  Mathematics,  Campus  Box 
170,  P.O.  Box  173364,  Denver,  CO  80217-3364.  bennethum@math.cudenver.edu.  To  whom 
correspondence  should  be  addressed. 

'Center  for  Applied  Math,  Math  Sciences  Building,  Purdue  University,  W.  Lafayette,  IN 
47907. 


1 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

08  AUG  2001  2' REPORT  TYPE 

3.  DATES  COVERED 

00-00-2001  to  00-00-2001 

4.  TITLE  AND  SUBTITLE 

Multicomponent,  Multiphase  Thermodynamics  of  Swelling  Porous  Media 
with  Electroquasistatics:  I.  Macroscale  Field  Equations 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Colorado  at  Denver, Center  for  Computational 

Mathematics, PO  Box  173364, Denver, CO, 80217-3364 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

see  report 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

18.  NUMBER  19a.  NAME  OF 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE 

unclassified  unclassified  unclassified 

32 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Swelling  Porous  Media  with  Electroquasistatics 


2 


multi-component  material.  We  combine  mixture  theory  and  averaging  to  obtain 
macroscopic  field  equations  for  a  deformable  porous  medium  in  which  species 
may  induce  different  electric  fields  but  forces  on  the  species  are  determined  by 
the  net  electric  field  of  the  mixture. 

The  constituent  forms  of  Maxwell’s  equations  were  first  developed  by  Kel¬ 
ly  [24],  although  he  did  not  consider  deformable  media.  A  binary  mixture  of 
electro-magnetic  fluids  has  been  used  to  model  plasma  with  single  electric  and 
magnetic  fields,  see  e.g.  Kulsrud  [25].  Benach  and  Muller  [2]  applied  single 
electric  and  magnetic  fields  to  model  a  mixture  of  dielectric  fluids.  The  idea 
of  using  multiple  electric  fields  has  been  used  historically  in  modeling  a  binary 
mixture  of  superconducting  fluids  [9].  Eringen  [14]  developed  a  mixture  theo¬ 
retic  approach  for  a  multiple  constituent  electro-magnetic  deformable  medium 
with  multiple  electro-magnetic  fields  and  applied  the  model  to  combinations  of 
conducting/superconducting  fluids  and  elastic  solids.  This  article  differs  from 
[14]  in  that  multiple  phases  are  considered,  a  different  philosophy  is  taken  in 
the  electric  fields  of  constituents  accounting  for  energy  and  momentum  balance, 
and  the  mixture  may  swell.  Mixture  theory  alone  does  not  allow  for  this  de¬ 
velopment.  An  upscaling  approach  is  required  to  map  microscopic  information 
into  well-defined  macroscopic  variables. 

On  the  microscale  one  can  distinguish  between  phases  or  materials,  while 
the  macroscale  is  defined  to  be  the  scale  at  which  the  material  appears  to  be 
homogeneous.  There  are  several  upscaling  techniques  to  choose  between:  ho¬ 
mogenization  (matched  asymptotics)  [5,  12,  35],  volume  averaging  in  the  sense  of 
Whitaker  [31,  33,  34,  41,  42],  spectral  integral  methods  [17],  generalized  Taylor- 
Aris  methods  [7],  and  hybrid  mixture  theory  [1,  3,  21,  22].  All  but  the  latter 
of  these  methods  upscale  field  equations  and  constitutive  equations  from  the 
microscale  to  the  macroscale.  The  hybrid  mixture  theoretic  approach  that  we 
adapt  does  not  upscale  constitutive  relations. 

We  present  constituent,  or  species,  electro-quasistatic  equations  for  a  mix¬ 
ture,  give  the  relationship  between  the  species  and  single-phase  properties,  and 
then  volume  average  these  governing  equations  for  a  multicomponent,  multi¬ 
phase  medium  to  obtain  the  governing  equations  for  the  medium  at  the  macroscale. 
At  this  scale,  the  medium  is  viewed  as  a  continuum  where  thermodynamic 
properties  for  each  constituent  of  each  phase  exist  spatially  everywhere.  The 
relationships  between  the  microscale  variables  and  the  macroscale  variables  are 
explicitly  given.  For  ease  of  exposition  we  assume  that  interfacial  properties 
such  as  excess  mass  density,  free  charge  on  interfaces,  and  interface  currents, 
are  negligible;  although  the  present  theory  can  be  extended  to  incorporate  these 
effects  [19].  The  medium  is  referred  to  as  multi-phase,  but  it  is  understood  that 
this  includes  composites  (e.g.  a  medium  composed  of  two  solid  materials),  or  a 
porous  medium,  where  the  multiphases  may  be  gas,  immiscible  liquids,  and/or 
solid. 

Maxwell’s  equations  involve  species  electric  fields,  but  the  macroscale  mo¬ 
mentum  and  energy  equations  for  species  are  based  on  the  philosophy  that  the 
force  or  work  induced  by  the  electric  field  acting  on  a  species  is  generated  by 
the  total  electric  field.  This  differs  from  what  has  been  done  in  the  past  [14,  24]. 


Swelling  Porous  Media  with  Electroquasistatics 


3 


The  species  electric  fields  are  defined  through  Gauss’  Law  since  the  charge  and 
polarization  density  are  well  defined  for  a  species.  The  species  electric  fields 
then  manifest  themselves  in  the  macroscale  equations  through  terms  involving 
gradients  or  time-rate-of-change  of  the  volume  fraction. 

The  resulting  field  equations  may  be  used  to  obtain  restrictions  on  admissible 
constitutive  relations  by  applying  the  methods  of  either  extended  thermodynam¬ 
ics  [30],  or  rational  thermodynamics  [13,  10].  The  need  for  these  equations,  and 
the  explicit  relationships  between  the  macroscopic  and  microscopic  variables  is 
demonstrated  by  the  increasingly  complex  porous  media  being  studied.  Recent 
examples  include  Huyghe  and  Janssen  [23]  and  Gu  et  al.  [20]  who  both  use  a 
simplified  version  of  the  equations  developed  herein  to  model  an  incompressible 
porous  medium  composed  of  an  electrically  charged  solid  phase  saturated  with 
an  ionic  fluid. 


2  Microscale  Equations 

In  this  section  we  present  the  governing  field  equations  at  the  microscale.  The 
equations  include  mass  balance,  conservation  of  linear  and  angular  momentum, 
conservation  of  energy,  entropy  balance,  and  Maxwell’s  equations.  The  complete 
form  of  these  equations  for  a  bulk  phase  (single  phase  with  no  species)  are  given 
by  Eringen  [15]  and  Tiersten  [38],  and  here  we  follow  their  formulation  and 
notation.  The  derivation  of  these  equations  follows  from  the  particle  level  in 
the  spirit  of  Lorentz  [27]. 

It  is  assumed  that  the  dominant  field  source  is  the  electric  field  and  that 
velocities  are  small  compared  to  the  speed  of  light  (non-relativistic).  Follow¬ 
ing  Melcher  [28],  the  electroquasistatic  system  of  equations  are  obtained  by 
non-dimensionalizing  Maxwell’s  equations,  expanding  each  variable  in  a  Tay¬ 
lor  series  about  the  variable  representing  the  ratio  of  the  electromagnetic  wave 
transit  time  to  the  characteristic  time  of  the  problem,  and  taking  the  zero-order 
equations  as  the  quasi-static  formulation.  The  difference  between  the  electro¬ 
quasistatic  formulation  and  the  magneto-quasistatic  formulation  arises  in  the 
choice  of  the  normalizing  parameters:  in  the  former  a  reference  electric  field  is 
used,  and  in  the  latter  a  reference  magnetic  field  is  used. 

It  is  necessary  to  postulate  governing  field  equations  which  hold  for  each 
specie.  Following  Truesdell  [40],  we  adhere  to  the  following  principles:  (1)  all 
properties  of  the  mixture  must  be  mathematical  consequences  of  properties  of 
the  constituents;  (2)  to  describe  the  motion  of  a  constituent,  we  may  perceive 
it  as  being  isolated  from  the  rest  of  the  mixture,  provided  we  allow  properly  for 
the  actions  of  the  other  constituents;  (3)  the  motion  of  the  mixture  is  governed 
by  the  same  equations  as  is  a  single  body.  Further,  due  to  the  incorporation 
of  electric  fields,  some  assumptions  must  be  made  regarding  the  form  of  the 
momentum  and  energy  balance  laws.  It  is  assumed  that  (4)  the  “primitive”  form 
of  the  balance  laws  is  the  one  incorporating  the  work  and  force,  and  not,  e.g.  the 
electro-stress  tensors.  This  is  in  contrast  to  [8,  26]  in  which  the  electro-magnetic 
stresses  are  considered  the  primitive  quantity.  The  forms  are  equivalent  up  to  the 
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classical  Maxwell’s  equations,  but  are  not  equivalent  in  the  mixture  formulation. 
Hence,  until  experiments  prove  otherwise,  the  choice  is  purely  philosophical. 
The  last  principle  is  that  (5)  within  the  electrical  work  and  force  terms  in  the 
energy  and  momentum  balance  equations  respectively,  it  is  the  total  electric 
field  which  acts  on  the  species,  and  not  just  the  electric  field  associated  with 
the  species.  This  last  assumption  guarantees  that  principles  (1),  (2),  and  (3) 
are  not  violated. 

Although  the  form  of  the  equations  for  the  constituents  is  natural  for  some 
balance  laws  (conservation  of  mass,  conservation  of  charge  density),  it  leads  to 
un-intuitive  variables  in  other  equations,  e.g.  the  partial  stress  tensor.  However, 
we  do  not  want  to  mix  bulk  phase  equations  with  constitutive  equations,  and 
hence  formulate  constitutive  equations  for  all  field  equations.  To  account  for 
the  effects  of  the  other  constituents,  exchange  terms,  denoted  by  variables  with 
a  carrot  •,  are  introduced. 

The  species’  equations  can  then  be  summed  to  obtain  the  governing  equation 
for  the  phase,  and  the  relationship  between  the  species  properties  and  their  bulk 
phase  counterparts  are  then  obtained.  The  details  regarding  this  procedure  for 
the  field  equations  with  no  electric  or  magnetic  effects  can  be  found  in  Trues- 
dell  and  Toupin  [39]  or  Bowen  [6] .  For  the  parallel  development  of  Maxwell’s 
equations,  see  Kelly  [24]. 

We  now  present  the  governing  equations  for  a  mixture  in  a  single  phase 
at  the  microscale  assuming  electroquasistatics.  The  relationships  between  the 
species’  variables  and  the  bulk  phase  variables  are  obtained  by  summing  the 
field  equation  over  all  constituents  and  relating  the  variables  so  as  to  obtain  the 
field  equation  at  the  bulk  scale.  These  relationships  are  given  in  Appendix  B. 

Conservation  of  Mass 


D\p>) 

Dt 


+  (V  •  v-’)  =  fPr’ 


where  is  the  material  time  derivative  given  by 


D-? 

15i 


■  V 


(1) 

(2) 


and  where  r°  is  the  rate  of  mass  transfer  to  species  j  from  other  species  due 
to  chemical  reactions.  Summing  over  constituents  we  obtain  the  bulk  phase 
equation 

^t+p(V-v)  =  0  (3) 

with  the  restriction 


N 

pip1  =  o. 

3=1 


(4) 


This  restriction  merely  states  that  within  an  isolated  system  consisting  of  a 
single  phase  there  is  no  net  loss  of  mass. 


Swelling  Porous  Media  with  Electroquasistatics 


5 


Gauss’  Law 

For  this  equation  we  introduce  the  electric  field  for  species  j,  E3 .  Gauss’ 
law  for  constituent  j  can  then  be  expressed  as 

V  ■  D1  -  ql  =  7  (5) 

where  D3  is  the  displacement  vector  for  the  jth  component  only,  q3e  is  the  charge 

density  of  constituent  j,  and  d!  is  the  excess  charge  density  due  to  the  presence 
of  other  species.  The  displacement  vector  for  constituent  j  is  defined  in  terms 
of  the  electric  field  and  polarization  density  as: 

DJ=e0Ej  +  Pj,  (6) 

where  Eq  is  the  permittivity  in  a  vacuum.  In  MKS  units,  the  permittivity  has 
the  value  of  8.854  x  10-12  Farads  per  meter  [28]. 

The  species’  electric  field,  E3 ,  must  satisfy  the  restriction  that  the  total 

N 

electric  field,  E ,  is  the  sum  of  the  species’  electric  field,  i.e.  J2E3  =  E.  If 

o= 1 

there  is  no  external  electric  held  then  E  is  just  the  electric  held  generated  by 
all  the  species.  There  is  no  unique  way  of  incorporating  the  externally  applied 
electric  held  into  these  equations.  One  could  treat  the  external  electric  held  as 
the  presence  of  another  species,  say  species  N,  in  which  this  “external”  species 
has  no  other  electro-thermodynamic  properties.  Then  species  N  would  have  no 
charge  density,  no  associated  polarization,  etc.  and  these  terms  would  be  set  to 
zero  in  the  governing  equations.  Alternately,  a  portion  of  the  external  held  could 
be  assigned  to  each  E3  so  that  E3  is  the  sum  of  the  electric  held  generated  by 
species  j  and  a  weight,  w3 ,  times  the  external  electric  held.  The  weights  must 
sum  to  one,  and  may,  for  example,  be  proportional  to  the  amount  of  charge, 
or  maybe  the  mass  fraction  of  species  j.  A  third  choice  might  be  to  define 
the  species  electric  held  so  that  the  exchange  term  7  is  zero.  This  simplifies 
computations  down  the  road  although  its  physical  interpretation  is  not  clear. 
These  choices  have  no  affect  on  the  following  derivation,  as  long  as  E3  sum 
to  the  total  electric  held.  However,  in  formulating  a  mathematical  model  it  is 
necessary  to  choose  a  particular  definition  and  remain  consistent. 

Summing  over  constituents  gives  the  bulk  phase  Gauss  law: 

V-D-qe  =  0  (7) 

1 v 

where  D  =  D3  and  with  the  restriction 

j= i 


N 

E^  =  °-  (g) 

3= 1 

which  states  that  the  net  effect  of  excess  charge  densities  produced  by  the  species 
acting  on  each  other  must  sum  to  zero. 
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Faraday’s  Law 

Faraday’s  law  for  the  quasi-static  case  for  species  j  is  given  by 

V  x  E3  =  o3  or  E3k  l  -  Ejk  =  £ikrno3m,  (9) 

where  the  second  form  of  the  equation  is  in  indicial  notation  with  repeated 
indices  denoting  summation,  a  comma  denotes  differentiation,  and  £ikm  is  the 
permutation  tensor.  Here  <rJ  incorporates  the  effect  of  the  electric  fields,  E\ 

Summing  over  species  gives  the  bulk  phase  Faraday’s  law, 

Vx£  =  0,  (10) 


N 

where  E  =  E3  is  the  total  electric  field.  This  relationship  is  assumed 
j— i 

throughout.  The  required  restriction  obtained  from  this  summation  is  that 

N 

E  = 0  (n) 

3=  1 


Ampere’s  Law 

Ampere’s  law  for  the  quasi-static  case  for  species  j  is  given  by 

J3  =  +  V  x  H3  -  V  x  ( P 3  x  V3)  +  h3  (12) 

— j 

where  H 3  is  the  magnetic  field  intensity  and  h  accounts  for  the  effects  of  other 
constituents. 

Summing  over  species  gives  the  bulk  phase  version  of  this  law, 

J=-^  +  VxcH-Vx(Pxd)  (13) 

where 


N 

E^  =  °-  (14) 

3= 1 

Conservation  of  Electric  Charge 

This  equation  can  be  derived  by  taking  the  divergence  of  Ampere’s  law  and 
the  time  derivative  of  Gauss’  Law  and  summing  the  results.  Of  all  Maxwell’s 
equations,  this  equation  is  the  most  accepted  in  mixture  form  since  each  com¬ 
ponent  has  a  well  defined  physical  interpretation.  The  conservation  of  electric 
charge  for  species  j  is 

V  •  (J3  +  q3ev3)  +  ^r=QJ  +  pW, 


(15) 
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where  v3  is  the  velocity  of  species  j  relative  to  a  fixed  frame  of  reference,  O’ 3 
is  the  free  current  density  of  constituent  j  measured  relative  to  species  j,  z° 
is  the  charge  per  unit  mass  and  q3  is  the  rate  of  gain  of  charge  density  due  to 
the  presence  of  other  constituents  but  not  due  to  chemical  reactions.  The  free 
current  density,  denoted  by  J~3 ,  is  related  to  the  free  current  density  relative  to 
a  fixed  (Eulerian)  frame  of  reference,  J3  ,  by 

J!  J‘  q!vJ.  (16) 

Another  form  of  this  equation  is  obtained  by  subtracting  out  the  conservation 
of  mass  to  reduce  redundancy,  and  this  yields 

■  TP  zJ 

P’^  +  V-J1  =?■  (17) 

Summing  over  constituents  yields 

V-(J  +  qev)  +  ^=  0,  (18) 

where 

N 

y]  [q3  +  p3 z3r3]  =  0.  (19) 

j= 1 

This  equation  states  that  the  net  gain  or  loss  of  total  charge  of  constituent  j 
due  to  ion  transfer  or  chemical  reactions  is  zero  in  an  isolated  system. 

In  addition,  because  of  the  coupling  between  Gauss’  law,  Ampere’s  law,  and 
the  conservation  of  charge,  there  is  a  coupling  between  the  exchange  terms: 

V  •  tiJ  -  =  q3  +  q3er 3 . 


Linear  Momentum  Balance 
This  equation  is  given  by 

fp  -  V  ■  t?  -  pjg3  -  q>eE  -  P‘  ■  VE  =  pj7  (20) 

where  t3  is  the  partial  Cauchy  stress  tensor  and  g3  is  the  external  body  force 

'TO 

acting  on  constituent  j.  The  exchange  term,  i  ,  takes  into  account  all  gain  of 
momenta  due  to  the  presence  of  other  species  but  not  due  to  chemical  reactions. 
The  last  term  on  the  left-hand-side  is  usually  referred  to  as  the  Kelvin  force, 
and  the  adjacent  term  (on  the  LHS)  is  referred  to  as  the  Lorentz  force  [15].  Note 
that  unlike  [14,  24]  we  assume  it  is  the  total  electric  field  which  contributes  to 
these  forces. 

Summing  over  constituents  gives  the  conservation  of  linear  momentum  for 
the  entire  phase, 


p—  -  V  t-  pg  -  qeE  -  P  ■  VE  =  0. 


(21) 
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In  this  form  t  is  the  Cauchy  stress  tensor  and  is  related  to  the  partial  Cauchy 
stress  tensors,  tJ ,  in  the  same  way  as  the  purely  thermo-mechanical  mixture  the¬ 
ory  (no  electric  field) .  This  is  what  naturally  appears  in  the  energy  conservation 
equation  given  below,  i.e.  if  other  terms  are  incorporated  into  the  definition  of 
the  bulk  phase  stress  tensor,  then  these  additional  terms  must  be  subtracted  out 
where  the  stress  tensor  appears  in  the  conservation  of  energy  and  the  restriction 
must  be  modified.  The  specific  relationship  between  the  constitutive  variables 
and  the  bulk  phase  variables  are  given  in  Appendix  B.  The  restriction  for  the 
conservation  of  linear  momentum  is: 

JV 

[p^i  +  =  0  (22) 

3= i 

The  restriction  states  that  the  net  momentum  gained  between  species  due  to 
mass  transfer  or  mechanical  momentum  must  be  zero. 


Angular  Momentum  Balance 

In  indicial  notation,  the  angular  momentum  equation  is 

CkLmt^ki  T  ChlmPk-^l  =  ^rn ’  (23) 

where  rhJm  is  the  ?nth  component  of  the  net  rate  of  gain  of  angular  momentum 
due  to  the  presence  of  other  species.  The  negative  sign  on  the  right  hand  side  is 
in  keeping  with  the  convention  that  exchange  terms  (in  this  case  rnJ )  represents 
a  rate  of  gain  of  a  property. 

Summing  over  constituents  yields 

Cklmt'kl  T  Cklm-PkPl  —  Om  (24) 


with  the  restriction  that 

N 

pi  fri:l  =  0.  (25) 

3  =  1 

This  equation  states  that  the  net  internal  angular  momentum  for  this  single¬ 
phase  medium,  consisting  of  the  sum  of  angular  momenta  generated  by  species 
acting  on  each  other,  must  be  zero. 

Conservation  of  Energy 

The  species  energy  balance  is 


r - 

Dt 


-  tj  :  Vvj  -  V  •  q  pjhj  -  J  -  E 


dPj 


-  V  •  (vjPj 


dt 

E  =  p*  & 


(26) 


where  the  colon  indicates  a  tensor  dot  product  (e.g.  a:b  =  JA  .  atJ btj )  and  QJ  is 
the  rate  of  energy  gain  due  to  the  presence  of  other  constituents  but  not  due  to 
chemical  reactions  (mass  transfer)  or  momentum  transfer.  The  sum  of  the  last 
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three  terms  on  the  left-hand-side  is  the  electrical  energy  source.  As  in  the  linear 
momentum  balance  we  assume  that  electric  energy  is  due  to  the  total  electric 
field,  not  just  part  of  it.  These  three  terms  take  on  a  variety  of  forms  [15,  28] 
which  are  related  through  the  Maxwell  equations  (Gauss’  law,  Faraday’s  law, 
Ampere’s  Law,  and  Conservation  of  Electric  Charge).  The  form  presented  here 
is  the  one  derived  directly  from  microscale  electrical  forces  using  a  statistical 
averaging  approach  (see  Eringen  [15])  making  no  use  of  Maxwell’s  equations, 

since  otherwise  exchange  terms  (cf\  <xJ  ,  <p)  appear. 

Summing  over  constituents  yields 

De  BP 

p—  -t:Vv-V -q-ph-J-E-  —  -E-[V-  (t>P)]  -E  =  0  (27) 

with  the  restriction 


N 


E 


pi  Qi  +  fPvi  ■  i  +  (eJ  + 


.  tP  )pJ 


=  0. 


(28) 


Entropy  Balance 

The  entropy  balance  equation  is  one  way  of  representing  the  second  law  of 
thermodynamics: 


pi  - V  •  <\>i  —  pity  =  pipi  +  pi  A-7  (29) 

where  fjJ  is  the  rate  of  transfer  of  entropy  from  other  constituents.  Here  AJ  is 
the  rate  at  which  entropy  is  generated. 

Summing  over  constituents  yields 


with  the  restriction 


p°^  -V-(/>  —  pb=  pi  A 

(30) 

N 

| pifji  +  piyiri  =  0. 

(31) 

i= i 


It  is  postulated  that  the  total  rate  of  entropy  generation,  obtained  by  summing 
over  all  species,  is  non-negative.  All  relations  for  this  equation  remain  unaltered 
from  the  purely  thermo-mechanical  mixture  theory. 


3  Macroscale  Equations 

In  this  section  we  average  the  microscale  equations  and  obtain  equations  at  the 
macroscale.  At  the  macroscale,  each  thermodynamic  variable  is  defined  spatially 
everywhere,  so  that  if  there  are  3  phases,  with  each  phase  containing  N  con¬ 
stituents,  the  medium  is  viewed  as  2>N  overlaying  continua.  The  first  subsection 
discusses  the  notation  and  assumptions  required  to  upscale,  and  the  following 
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subsections  present  the  macroscale  field  equations  and  how  each  variable  relates 
to  the  microscale.  To  the  authors’  knowledge,  the  electroquasistatic  equations 
have  not  been  presented  before. 

3.1  Averaging  Procedure 

Consider  a  multi-constituent  multi-phase  medium  where  the  phase  is  denot¬ 
ed  by  Greek  letters  (a,  (3,  7),  and  the  constituent  or  species  is  denoted  by 
j,  j  =  1  For  ease  of  exposition  we  assume  that  interfaces  contain  no 

thermodynamic  properties.  Consequently  it  is  assumed  no  amount  of  mass, 
momentum,  energy,  charge,  current,  etc.  is  lost  when  being  transferred  between 
phases.  Interfacial  effects  can  be  included  by  pursuing  any  of  the  approaches  of 
[18,  19,  29],  however,  we  shall  omit  these  terms  to  keep  the  level  of  algebra  at  a 
minimum. 

The  governing  microscopic  equations  for  each  phase  were  given  in  Section  2, 
but  to  distinguish  between  phases  we  introduce  the  additional  Greek  superscript. 

Assuming  no  surface  discontinuities,  the  constituent,  microscopic  field  e- 
quations  of  mass,  charge,  linear  momentum,  angular  momentum,  energy,  and 
entropy  can  be  expressed  for  a  given  phase,  a,  as  (following  the  notation  of 
Eringen,  [13]): 

+  V  •  {p>vjft)  V  •  ij  ft  f 3  -  F3  =  ft  G3  +  pi  ft  (32) 

where  ip3  is  the  mass-average  (over  the  phase)  thermodynamic  property  of  con¬ 
stituent  j ,  v3  is  the  mass-average  velocity  vector,  ft  is  the  mass  density,  i3  is  the 
flux  vector,  f3  is  the  body  source,  F3  is  the  supply  due  to  electrical  effects,  G3  is 
the  net  production,  and  ip  represents  the  influx  of  ip  from  all  other  constituents 
(e.g.  due  to  chemical  reactions).  If  there  is  only  one  constituent,  %p  is  zero. 
For  each  of  the  respective  equations,  the  quantities  given  in  Table  1  are  used. 


Table  1:  Quantities  for  Equation  (32) 


Quantity 

%p 

i 

f 

F 

i3 

G 

Mass 

1 

0 

0 

0 

r 

0 

Charge 

z 

j 

0 

0 

q  +  pzr 

0 

Linear  Momentum 

V 

t 

g  +  9i 

F 

i  +  rv 

0 

Angular  Momentum 

r  x  v 

r  x  t 

r  x  g 

r  x  Fe  +  Ce 

r  x  (i  +  rv) 

0 

Energy 

e+\y2 

tv  +  q 

g  ■  v  +  h 

we 

Q  +  i-  v  +  r(e  +  \ v 2) 

0 

Entropy 

T) 

<t> 

b 

0 

rj  +rv 

A 

For  conciseness  we  also  introduce  the  force,  couple,  and  work  due  to  the  electric 
field,  respectively: 

F3  = 

e 


q3eE  +  P3  VE 


(33) 
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=  -V  •  ( DjE  -  \e0Ej  ■  El)  +7  E  +  \e0  (VE  ■  Ej  -VEj  ■  E\ 34) 


Cl  =  P-i  x  E 


(35) 


wi 


.I  E  —  E  V  (v'E  ■  P‘) 

^  (~e0EJ  ■  E)  -  V  •  (E  x  W)  +  V  •  (PV  •  E) 


1  (dEP 

+K  'E~2£°\W 


-f  - 


(36) 


(37) 


(see  also  equations  (20,23,26)).  Definitions  of  all  terms  are  given  in  Appendix 
A.  The  forms  of  FJ  and  IT)?  are  a  consequence  of  which  form  one  assumes  to 
be  the  most  primitive.  If  one  assumes  that  the  second  forms  of  FJe  and  Wj, 
(34,  37),  are  the  most  primitive  then  miscellaneous  terms  would  appear  in  the 
first  forms.  The  two  forms  are  needed  in  order  to  obtain  conditions  which 
correspond  directly  with  jump  boundary  conditions  across  an  interface.  The 
averaging  procedure  is  based  on  ideas  laid  down  in  [16,  41,  42,  37].  Several 
methods  are  available,  but  we  choose  the  computationally  simplest.  Equations 
are  averaged  by  weighted  integration  using  the  indicator  function  of  the  a-phase. 
To  avoid  the  mathematical  difficulties  of,  for  example,  defining  a  derivative  of 
the  averaged  quantities  resulting  from  using  such  a  weighting  function,  one  must 
treat  the  averaged  quantity  as  a  distribution  [36,  32]. 

It  should  be  noted  that  using  this  simple  weight  function  may  mean  that 
the  averaged  value  may  not  represent  the  actual  values  being  measured.  To 
account  for  the  measuring  technique,  one  needs  to  choose  a  weight  function 
which  represents  the  instrument  used  to  measure  the  physical  properties  [11]. 
Extensions  of  the  presented  theory  to  such  cases  are  straight  forward. 

Let  SV  be  a  volume,  SVa  the  portion  of  SV  in  the  a-phase,  and  SAap  the 
portion  of  the  interface  within  SV.  It  is  assumed  that  SVa  and  5A[S  are  isolated 
simply  connected  regions.  If  the  magnitude  of  SV  is  denoted  by  |<5Vj  then  the 
volume  fraction  can  be  expressed  as 


£a 

ix  t)  -  l^»l 

\6V\ 

(38) 

so  that 

a 

(39) 

The  indicator  function  is 

7 oc{r,t)  =  < 

r  i 

l  0 

if  r  e  SVa 

if  r  €  SV/3,  (3  yf  a. 

We  may  write  r  =  x  +  £  where  x  is  the  macroscale  spatial  variable,  and 
£  varies  over  SV.  To  obtain  the  macroscale  equations  formally,  one  multiplies 
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equation  (2)  by  7a(r),  integrates  over  SV  with  respect  to  £  (x  is  held  fixed)  and 
divides  by  \5V\ .  Then  in  order  to  obtain  equations  of  the  forms  which  mirror 
the  microscale  equations,  the  following  theorem  is  applied  to  interchange  the 
order  of  differentiation  and  integration. 

Theorem  1  If  wap  is  the  microscopic  velocity  of  interface  a(3  and  na  is  the 
outward  unit  normal  vector  of  5Va  indicating  the  integrand  should  he  evaluated 
in  the  limit  as  the  a/3 -interface  is  approached  from  the  a-side  then 

m  =  i[w\SwM&. 

'  n*da(0  (40) 

/J#a  1  V]  JSA“0 


+^w\Ljn-da(()-  (41) 

After  averaging  equation  (32),  the  system  is  considered  to  be  a  mixture  so 
that  each  component  in  each  phase  and  each  bulk  phase  now  have  thermo¬ 
dynamic  properties  existing  at  each  point  within  the  macroscopic  body.  The 
macroscopic  definition  of  each  field  variable  in  terms  of  its  microscopic  counter¬ 
part,  making  no  small  perturbation  assumptions,  is  given  in  Appendix  C. 

For  more  details  regarding  this  procedure  for  the  field  equations  with  no 
electric  or  magnetic  effects  see  for  example  Whitaker  [41],  Slattery  [37]  and 
Plumb  and  Whitaker  [31].  For  the  parallel  development  of  Maxwell’s  equations 
with  constitutive  assumptions,  see  Rio  and  Whitaker  [33,  34]. 
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where  is  the  velocity  of  species  j  at  interface  a/3. 

The  bulk  phase  counterpart  is  given  by 

jjafFa  a\ 

- — - +  e  p  (V  •  v  )  =  2^  £  P  ep.  (45) 

The  net  gain  of  mass  of  the  bulk  phase  must  be  zero,  implying  that: 

N 

pajr°‘:i  =0  V  a.  (46) 

t=i 

Further,  since  the  interface  is  assumed  to  be  massless,  we  have  the  restriction: 

eapa +  e?pP*e%  =0,  j=  1, ...,  N.  (47) 

Using  (45)  we  can  re-write  (42)  as 


_<A  cx. 
£  p 


DaC° 

Dt 


+  V  •  (e  p  ]v  ■”  )  =  2^  £  P  J(e0  -e0)+£  p 


(48) 


Pt^OL 


Gauss’  Law 

The  macroscopic  form  of  Gauss’  law  for  constituent  j  in  phase  a  is 

V  •  (eaD“0  -  eaq^  =  £aTj  +  Y  ead*p  (49) 


where  D j  is  the  volume  average  of  D° .  Here  dp  represents  the  effect  con- 
stituent  j  in  phase  (3  has  on  the  charge  of  the  same  constituent  in  phase  a: 

(50) 

Other  relations  are  given  in  Appendix  C. 

Summing  over  j  yields  the  bulk  phase  form: 

V-(saDa)-eaq:=Y£aT0-  (51) 

0^a 


The  restrictions  include 


N 

Y^  =0  Va  (52) 

3=i 

eaT0  +  e0da  =0,  j  =  1, ...,  N.  (53) 

The  first  restriction  states  that  the  net  effect  constituents  have  on  each  other 
within  phase  a  must  be  zero,  and  the  second  states  that  the  interfacial  displace¬ 
ment  is  zero,  which  is  a  consequence  of  assuming  that  the  interface  contains  no 
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electric/thermodynamical  properties.  Recall  that  if  there  is  no  surface  charge 
or  polarization  the  jump  condition  between  two  materials  states  that  the  jump 
in  the  normal  component  of  the  displacement  vector  must  be  zero  [15]: 

n  ■  [[£>]]  =  0.  (54) 

Integrating  this  equation  over  6Aap  shows  that  this  theory  is  consistent  with 
the  classical  formulation  for  the  jump  condition,  equation  (53). 

Faraday’s  Law 

The  macroscopic  form  of  Faraday’s  law  for  constituent  j  in  phase  a  is 

V  x  (eaEai)  =  £aoa>  +  Y  eaaa0j  (55) 

/3^a 


or  in  indicial  notation 


(eaE*i),l-(eaE^),k=elkmead0ri 


m 


(56) 


where 


<?  a  =  ~ 


\SVa\  J§Aa 


na  x  E J  da. 


(57) 


which  represents  the  effect  phase  (3  has  on  the  electric  field  of  phase  a.  Upon 
summing,  the  bulk  phase  form  of  Faraday’s  law  becomes: 


V  x  ( eaEa )  =  Y  £aV0 

0^a 


(58) 


with  restrictions 

N 

Yj  & 0,3  =  0  V  a 

(59) 

3  =  1 

eaa^  +  £0$%  =0  j  =  1, ...,  N. 

(60) 

Equation  (60)  corresponds  precisely  with  the  classical  jump  condition  across 
interfaces  [15] 


n  x  [[£]]  =  0,  (61) 

where  the  double  square  brackets  indicate  the  difference  of  the  quantity  evalu¬ 
ated  on  either  side  of  the  interval,  if  (61)  is  integrated  over  the  surface  SAap. 
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Ampere’s  Law  The  macroscopic  form  of  Ampere’s  law  for  the  quasi-static  case 
for  species  j  of  phase  a  is  given  by 


eaJai  =  - 


d(eaDa J) 
9t 


V  x  (eaHa*) 


V  x  (eaPaj  x  va*)  +  53  ^7  +  (62) 


where 


J5Aa 


D3w3afi  ■  na  -  ( HJ  -  PJ  x  vJ)  x  n° 


da 


(63) 


represents  the  effect  of  phase  j3  on  phase  a.  Summing  over  species  gives  the 
bulk  phase  version  of  this  law, 


eaJa  =  - 


d(eaDa) 

dt 


V  x  eaHa  -  V  x  ( saPa  x  va)  +  53  £ahp 


where  the  restrictions  are 

N 

Y  h  1  =  0  Va 

3  =  1 

£ahap  +  =  0  j  =  l,...,N. 

The  jump  condition  associated  with  Ampere’s  law  is  given  by  [15]: 
n  x  [[H  —  w  x  D]\  =  0 


(64) 


(65) 

(66) 


(67) 


where  the  polarization  at  the  surface  has  been  neglected.  Equation  (63)  does 
not  correspond  directly  with  the  jump  condition  given  in  [15]  since  (67)  is  de¬ 
rived  directly  from  the  global  surface  area  form  of  Ampere’s  law,  and  (63)  is 
derived  from  the  bulk-phase  form  of  the  law.  This  causes  a  loss  of  information; 
specifically,  (63)  does  not  include  a  term  corresponding  to  wJ  (n  ■  DJ). 

Conservation  of  Electric  Charge 

The  conservation  of  charge  equation  at  the  macroscale  becomes 


V  •  (eaJa>  +  £aq°*va!)  +  =  eaT3  +  £apa3z^raj 


+  53  eapa>(Zap1 +za%’) 


(68) 


where 


20  Pa3\SVa\JsAa 


■  na  da 


Pa°  \&Va\  JSAa 


(wa0  ~  v’’)  '  da  (69) 
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which  represents  the  rate  of  exchange  of  charge  of  constituent  j  from  phase 
(3  to  phase  a  not  clue  to  mass  exchange.  Using  the  continuity  equation,  (42), 
equation  (68)  may  be  re-written  as: 


eapa’  °a^3  +  V  •  {eaJa*)  =  +  X]  £C*Pajz7 


Dt 


Summing  over  constituents  yields 


07 ^ 


v  •  (eaJa  +  eaq“va)  +  ^(e“(£)  =  £  sapa(Za0  +  zae%), 

0+01 


where  the  following  restrictions  apply 

N 

ft*3  +  Paj z^r011]  =0  Va 
f=i 


eapa’(Zj  +  z^ep)  +  e0^(zij  +  z03e03)  =  0  j  =  1,  ...,1V.  (73) 


(70) 


(71) 


(72) 


Equation  (73)  corresponds  precisely  with  the  classical  jump  condition  across  a 
discontinuous  interface  [15]: 


n  ■  [[J  +  qe(v  -  w)]}  =  0. 


(74) 


Linear  Momentum  Balance 

The  macroscale  linear  momentum  equation  is  given  by 


eapa 


Da j  Vai 

Dt 


--£qEt  •  Eaj'Ve° 
2 


V  •  (eaia0  -  eapa3  {g01*  +  g “*)  -  £aq^ET 
.  e“p“i  •  VEt  =  £apa>TJ  +  £<xpaj  t7 

0+ot 


(75) 


where  Et  is  the  upscaled  (total)  electric  field  and  is  related  to  the  species 

N 

electric  field  by  Et  =  and 

1=1  a 


Taj  1  [ 

0  ~  Paj\6Va\  JSAa 


(■ tj)T  +  ED1  -  -£0E3  ■  El  +  p3v3(w3a0  -  v3) 


na  da 


Va j  C  1 

- n — r  /  p3  (w3  R  —  v3)  ■  na  da -\ - d3  Et 

Pa3\SVa\JsAa0  K  a 0  Pa*  0 


represents  the  effect  constituent  j  of  phase  (3  has  on  the  rate  of  change  of  me¬ 
chanical  momentum  of  the  same  constituent  in  phase  a.  The  relation  between 
the  other  macroscale  variables  and  their  microscale  counterparts  are  given  in 
Appendix  C. 


(76) 


Swelling  Porous  Media  with  Electroquasistatics 


17 


There  are  a  couple  of  additional  forces  in  this  form  compared  to  the  classical 
form  due  to  the  introduction  of  the  species  electric  field,  Eaj ,  and  the  volume 
fraction  ea.  The  first  we  consider  an  additional  (internal)  body  force,  gj.  It 
is  due  to  the  difference  between  E  and  E3  (see  Appendix  C),  and  so  if  the 
source  of  the  electric  held  is  dominated  by  one  charged  species,  then  this  term 
is  negligible.  The  other  additional  term  involves  the  gradient  of  the  volume 
fraction.  If  the  medium  is  homogeneous  in  volume  fraction  ( ea  is  constant), 
then  this  term  is  zero.  The  form  of  this  equation  is  motivated  by  the  form  of 

a  j 

T p  ,  which  we  choose  to  be  consistent  with  the  jump  condition  associated  with 
momentum  balance. 

Summing  over  j  yields 


zapaD-j^-  -  V  •  (e“t“)  -  +  9?)  -  eaQeET  -  eaPa  ■  VET 


-e0ET  •  EaVea  =  J2  zapaTp,  (77) 


with  restrictions 

N 


E  paj  (i  J  +  r aj  vaj )  =  0  V  a 
j= i 

•  a  1  '-'a. 


(78) 


£>^(T:J+ff^)+EV(^+e^ft)=  0  j  =  1, ..., N.  (79) 


P  'rc/3 

The  jump  condition  across  a  discontinuous  interface  is  [15]: 

n  ■  [[pv(v  —  w)  —  t  —  tE}}  =  0  (80) 

where  w  is  the  speed  of  the  discontinuity  and  for  the  electroquasistatic  case 
considered  here, 


Fe  =  V-tE  =  V-  (DE  ^eoE-EI).  (81) 

Equation  (80)  can  be  shown  to  correspond  directly  with  (79)  using  (53). 
Angular  Momentum  Balance 

The  macroscale  form  of  the  conservation  of  angular  momentum  equation  in 
indicial  notation  is 


saekimQ  -  £aeklmP^(ET)i  =  ea pajfh°ri 


+  E  £apaj 

p^ot 


+  eapaiM° 


(82) 


where  the  definitions  of  all  variables  are  given  in  Appendix  C.  Thus  the  macroscale 
equation  is  composed  of  two  parts.  The  equation  consisting  of  all  but  the  term 
involving  Maj  is  the  macroscale  form  of  the  microscale  equation.  The  term 
Maj  contains  all  terms  arising  from  the  the  microscale  conservation  of  linear 
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momentum  crossed  with  £  and  volume  averaged.  These  additional  terms  for  the 
purely  mechanical  case  have  been  derived  before  ([16,  22]).  The  term  in  Maj 
involving  the  time  rate  of  change  of  velocity  is  known  as  inertial  [16]  or  local 
[22]  spin.  The  term  in  Maj  involving  t3  has  been  referred  to  as  the  first  sur¬ 
face  stress  moment  [16],  or  apparent  couple  stress  tensor  [22].  One  observation 
which  is  painfully  clear  is  that  even  in  a  medium  which  contains  no  electric  field 
or  charge,  the  mixture  of  two  media  which  have  symmetric  stress  tensors  may 
result  in  a  macroscopic  medium  which  has  a  non-symmetric  stress  tensor. 

Summing  over  j  yields 

-eaeklmth  ^  £aeklmP?(ET)i  =  £apa’M«  +  ^  eapa{rr%)m  (83) 

where  we  note  that  symmetry  of  the  stress  tensor  is  lost  due  to  possible  polar¬ 
ization  of  the  medium,  the  microscale  angular  momentum,  and  the  interaction 
of  the  phase  with  other  phases. 

The  restriction  which  must  hold  to  preserve  angular  momentum  within  a 
phase  is 

N 

0  V  a.  (84) 

j=i 

The  restriction  which  arises  from  the  assumption  of  no  interfacial  angular  mo¬ 
mentum  is 

£apa’my  +  =0  j  =  1, ...,  N.  (85) 


Conservation  of  Energy 


The  conservation  of  energy  equation  is  given  by 

eV 


T)Oij  ~otj 

-£TbVi)a>  -V  •  {eaqaj)-  ea paJ ha>  -  ea  Jaj  ■  ET 


Dt 
d(£aPai) 
dt 


1  Daj£a 

Et  —  V  •  (eava3Pa3)  •  Et  -  -£0 Et  ■  Ea* - 

v  )  2  Dt 

=  £apa>  Qa>  +  zapa°  Q°p  (86) 

0  ±OL 


where 

Q7  = 


pmti  L  >J + (t' + D‘E  -  ie°(E’ ' E)I) ' 


(3 

+p3(e3  +  ^v3  •  v3  +  ^ e0 E3  ■  E)(w3a/3  —  v J')]  •  na  da 


eai  -  \va>  ■  vai 


paj  \sva 


lSAa 


pj{wip  —  v3)  ■  na  da 


PaiWa\  JSAC 


O  L 


( tj)T  +  ED3  -  -£0(E3  ■  E)I  +  p3v3(w3a0  -  v 3) 


■  na  da 
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+^g  La  e [( “ £oEJE  +  £o{E1 ' E)I) ' ^  +  H>  xE]'na  da 

— Ot-i  ^ CXj 

-d/  Et  ■  va>  -  hp  ■  Et  (87) 


represents  the  effect  constituent  j  of  phase  j3  has  on  the  rate  of  change  of  energy 
of  the  same  constituent  in  phase  a  due  to  non-mass  transfer,  non-mechanical 
means.  The  relation  between  the  other  macroscale  variables  and  their  microscale 
counterparts  are  given  in  Appendix  C.  As  in  the  linear  momentum  equation, 
upscaling  produces  an  additional  term  involving  the  material  time  rate  of  change 
of  the  volume  fraction.  This  term  would  be  negligible  in  a  non-swelling  porous 
medium.  Additional  effects  of  fluctuations  in  the  species  electric  work  terms 
manifest  themselves  in  the  flux  term,  qaJ ,  and  the  external  source  term,  h°j . 

Summing  over  constituents  yields 


Daea  d(saPa 1 

-  eata  :  Vva  -  V  •  ( £°qa )  -  £apaha  -  £aJa  ■  Et  -  1  ’  ■  ET 


£  p 


Dt 


dt 

-V  •  ( £avaPa )  ■  Et  -  l-£oEa  •  Et  ^  eapaQa0  (88) 


with  the  restrictions 
N 


E 

3= 1 


pa:i  Qa:i  +  paH  j  ■  va>  +  pa ^  (  eaj  +  -(u“02 


=  0.  V  a 


(89) 


£apaiQ°p  +  £ap“JT/  •  vai  +  £apa>ep>  (  eaj  +  -(«“') 


£0P0jQa  +  £PpPjTa  ■  V +  £0p0je%  (  +  ~(v&)2 


=  0  j  =  l,...,N. 

(90) 


The  first  restriction  states  that  energy  must  be  conserved  within  a  phase,  and 
the  second  states  that  the  interface  can  hold  no  energy. 

We  wish  to  compare  (90)  with  the  jump  discontinuity  condition.  Various 
forms  exist  depending  on  what  is  considered  negligible  and  the  use  of  Maxwell’s 
equations.  Here  we  consider 

[[q  +  t  ■  v  +  (pe  +  ^pv2){w  —  v)  +  ^£$(E  ■  E)w  +  H  x  E  +  PE  •  «]]  ■  n  =  0.  (91) 
Equation  (90)  can  be  re-written  as: 


f  q 3  +  (tj+PjE)-  vJ  +  p?  {e?  +  ) 2 )  (wjafj  -  vj ) 

l  * 


\SV\  J8A, 


|£V|  IsAc 


-~£0{&  -E)w^  +  W  xE 


■na  da 


«  L 


+  (tJ  +  P-i E)  ■  iV  +  pJ(eJ  H — (v^)2)(wap  ~  v J) 
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--eo(E^E)w^- 


H3  x  E 


i13  da.  (92) 


Comparing  terms  we  see  that  the  condition  involving  exchange  of  energy 
across  the  interface,  (90),  conforms  exactly  with  the  classical  jump  condition 
(91). 

Entropy  Balance 

The  entropy  balance  for  constituent  j  in  phase  a  is  given  by 


Daj  naj 
£apa*  ' 


Dt 


-V-ea0“'-e°79a'&a*  =  ^£apa^+£ap0‘3rfl3+£°‘pa3ha 

(93) 


where  represents  the  rate  at  which  entropy  is  gained  from  constituent  j 
in  phase  (3.  Summing  over  constituents  we  get 

£apa^f-  -  V  •  (£“</>“)  -  eapaba  =  Y  £aPa®0  +  £“p“AQ,  (94) 


Dt 

with  restrictions 


N 


'Y ,  Paj  {vf*3  +  r“J  r]aj )  =  0  V  a. 

3=1 


(95) 


£ap“J(V  +  3SV0  +eV'(*a  +e%riY  =0  j=  1  (96) 


The  second  condition  is  exactly  analogous  to  the  classical  jump  boundary  con¬ 
dition  applied  at  an  interface.  This  equation  can  be  used  to  obtain  restrictions 
on  the  forms  of  constitutive  equations  using,  for  example,  the  Coleman  and  Noll 
method  (see  e.g.  [4,  13]).  At  this  point  we  make  no  assumptions  regarding  the 
form  of  the  entropy  flux  or  source.  It  is  now  postulated  that  the  rate  of  entropy 
production  is  non-negative  after  summing  over  both  species  and  phases.  Our 
statement  at  the  end  of  Section  2  need  not  hold  in  this  case,  as  it  is  the  rate  of 
entropy  production  of  the  universe  which  is  non-negative. 


4  Discussion 

We  have  provided  the  microscale  and  macroscale  balance  laws  for  a  multiphase, 
multi-species,  swelling  system  which  incorporates  electroquasistatics.  The  mi¬ 
croscale  derivation  follows  the  approach  of  Eringen  and  Maugin  [15],  while  the 
upscaling  to  the  macroscale  is  achieved  via  spatial  averaging.  The  species  elec¬ 
tric  field,  Eaj ,  is  carefully  defined  and  the  complete  relations  between  the  mi¬ 
croscopic  and  macroscopic  variables  and  the  species  and  bulk  phase  variables 
are  provided.  Further  comparisons  were  made  between  restrictions  required  to 
hold  at  interfaces  and  jump  conditions  classically  used  at  interfaces. 
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In  regards  to  the  conservation  of  momentum  and  energy,  there  are  several 
ways  in  which  one  can  group  together  terms  into  the  definitions  of  the  macro¬ 
scopic  variables.  In  this  paper,  we  choose  the  macroscopic  definitions  in  a  way 
which  results  in  consistency  between  restriction  conditions  at  interfaces  and 
classical  jump  conditions  used  at  interfaces.  But  as  a  consequence,  additional 
body  sources  appear  due  not  to  external  sources  but  due  to  species  electric  fields 
and  microscopic  fluctuations  of  the  work  terms.  Further,  there  are  additional 
terms  which  appear  due  to  gradients  in  volume  fractions  and  material  time  rate- 
of-change  of  the  volume  fractions  in  the  conservation  of  momentum  and  energy 
equation,  respectively.  This  is  not  due  to  the  primitive  choice  of  FJe  and  W£  (see 
equations  (33-35)).  Rather  it  is  a  consequence  of  the  choice  of  grouping  terms 
so  as  to  recover  classical  jump  conditions.  The  consequence  of  these  additional 
terms  needs  to  be  investigated.  This  is  partially  done  in  part  II  of  these  papers. 

The  field  equations  are  required  for  many  practical  applications  including, 
but  not  limited  to,  swelling  clays  (smectites),  biopolymers,  biological  mem¬ 
branes,  pulsed  electrophoresis,  and  chromatography.  In  part  II  of  these  papers 
we  exploit  the  entropy  inequality  with  independent  variables  consistent  with 
several  natural  systems  and  apply  the  resulting  constitutive  theory  near  equi¬ 
librium  for  two  problems:  electrolyte  transport  in  swelling  clays  and  pulsed 
electrophoresis. 
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Appendix  A.  Nomenclature 

Superscripts,  Subscripts,  and  Other  Notations 

■ai  jth  component  of  a-phase  on  mesoscale 
•a  a-phase  on  mesoscale 
~  fluctuation  from  averaged  quantity 
•  denotes  exchange  from  other  interface  or  phase 
difference  of  the  two  quantities,  i.e.  • k  —  -l 
•L.  microscopic  property  of  constituent  j  in  phase  [subscript]  (non-averaged) 

Latin  Symbols 

5Aap:  Portion  of  a/3-interface  in  representative  elementary  volume  (REV) 
baj,ba:  External  entropy  source  [J/(Kg-s-°K)] 

Cai :  Mass  fraction  of  j^1  component  [-] 

DQj,Da:  Electric  displacement  [C/m2] 

— Otj 

d  :  Net  effect  other  constituents  have  on  the  charge  of  constituent  j  within 
phase  a  in  Gauss’  Law  [C] 

— Otj 

dp  :  Net  effect  phase  (3  has  on  the  charge  of  constituent  j  in  phase  a  in  Gauss’ 
Law  [C] 
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eaj,ea:  energy  density  [J/Kg] 

' e g3:  Rate  of  mass  transfer  from  phase  [subscript]  to  phase  [superscript]  per 
unit  mass  density  [1/s] 

Eaj ,  E°:  Electric  Field  intensity  generated  by  (the  jth  constituent  of)  phase 

a  [V/m] 

gaj  ,ga:  External  supply  of  momentum  (gravity)  [m/s2] 

g“J  ,gj:  Internal  supply  of  momentum  due  to  fluctuations  in  Kelvin  and 

Lorentz  forces  [m/s2] 


,  ha:  External  supply  of  energy  [J /  (Kg-s)] 

Otj  ^ ot 

h  ,h  :  Free  current  density  induces  by  other  species  within  the  same  phase 
within  Ampere’s  law  [A/m2] 

-^Otj  ^ Ot 

hp  ,  hp:  Free  current  density  induces  by  other  phases  on  species  j  in  Ampere’s 
law  [A/m2] 

"7 Otj 

i  :  Rate  of  momentum  gain  due  to  interaction  with  other  species  within  the 
same  phase  per  unit  mass  density  [N/Kg] 

Haj  ,Ha :  Magnetic  field  intensity  [A/m] 

Otj 

i  :  Rate  of  momentum  gain  due  to  interaction  with  other  species  within  the 
same  phase  per  unit  mass  density  [N/Kg] 

jaj .  ju.  pree  current  density  in  a  fixed  frame  of  reference  [A/m2] 

Ja\Jn-  Free  current  density  in  the  material  frame  of  reference  [A/m2]  = 
[C/m2-s];  J'x‘  =  Jaj  -  qVva> 

rnaj :  Rate  of  angular  momentum  gain  due  to  interaction  with  other  species 
within  the  same  phase  per  unit  mass  density  [N-m/Kg] 

m'p  :  Rate  of  angular  momentum  gain  by  constituent  j  in  phase  a  due  to 
interaction  with  phase  f3  [N-m/Kg] 

Maj :  Rate  of  angular  momentum  gain  due  to  the  microscale  angular  momen¬ 
tum  terms  -  see  Appendix  C  [N-m/Kg] 

na:  Unit  normal  vector  pointing  out  of  a-phase  within  mesoscopic  REV  [-] 

pai  _  pa-  Polarization  density  averaged  over  a-phase  [C/m2] 

qeJ  Charge  density  averaged  over  a-phase  [C/m3] 
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rnaj :  Rate  of  gain  of  angular  momentum  of  consitutent  j  from  other  con¬ 

stituents  in  phase  a  [m2/s2] 

Partial  heat  flux  vector  for  the  component  of  phase  [J/(m2-s)] 
qa:  Heat  flux  vector  for  phase  a  [J/(m2-s)] 

:  Net  rate  of  charge  density  gain  by  species  j  due  to  interaction  with 

other  species  within  phase  a  (does  not  include  that  gained  to  due  mass 
transfer).  [C/m3-s] 

Qaj :  Rate  of  energy  gain  due  to  interaction  with  other  species  within  the 
same  phase  per  unit  mass  density  not  due  to  mass  or  momentum  transfer 
[J/(Kg-s)] 

Q^3  ,  Q\ }■  Energy  transfer  rate  from  phase  [subscript]  to  phase  [superscript] 
per  unit  mass  density  not  due  to  mass  or  momentum  transfer  [J/(Kg-s)] 

r :  Microscale  spatial  variable  [m] 

rQj :  Rate  of  mass  gain  due  to  interaction  with  other  species  within  the  same 
phase  per  unit  mass  density  [1/s] 

t :  Time  [s] 

T:  Temperature  [°K] 

taj:  Partial  stress  tensor  for  the  j ^  component  for  phase  [N/m2] 
ta:  Total  stress  tensor  for  the  phase  [N/m2] 

^  otj  ^  a 

T p  ,Tp:  Rate  of  momentum  transfer  through  mechanical  interactions  from 
phase  [subscript]  to  phase  [superscript]  per  unit  mass  density  [N/Kg] 

vaj ,  Va'.  Velocity  [m/s] 

5V:  Representative  elementary  volume  (REV) 

SVa:  Portion  of  a-phase  in  REV 

wai3:  Velocity  of  constituent  j  at  interface  between  phases  a  and  /3  [m/s] 
x :  Macroscale  spatial  variable  [m] 

zaj:  Charge  per  unit  mass  density  of  constituent  j  in  phase  a  [C/Kg],  qf3  = 

pQLj  gOt-j 
—  Otj  ^ oc 

Z p  ,Zp\  Rate  of  exchange  of  charge  of  constituent  j  from  phase  /3  to  phase 
a  per  unit  mass  [C/Kg-s] 
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7a :  Indicator  function  which  is  1  if  in  mesoscopic  region  a  and  zero  otherwise 
ea:  Volume  fraction  of  a-phase  in  mesoscale  REV  [-] 

£(P  Permittivity  in  a  vacuum.  In  MKS  units,  the  permittivity  has  the  value 
of  8.854  x  10-12  Farads  per  meter. 

Ckim-  Permutation  tensor 

AQj  ,  A“:  Entropy  production  per  unit  mass  density  [J/(Kg-s-°K)] 

Microscale  spatial  variable  which  varies  over  REV  for  fixed  x:  r  =  x  +  £ 

[m] 

T7aV?7Q:  Entropy  [J/(Kg-°K] 

rjaj :  Entropy  gain  due  to  interaction  with  other  species  within  the  same 
phase/interface  per  unit  mass  density  [J/(Kg-s-°K)] 

< f)aj :  Partial  entropy  flux  vector  for  the  j ^  component  for  phase  [J / (m2-s-°K)] 

< pa :  Total  entropy  flux  vector  for  the  phase  [J/(m2-s-°K)j 

,4>^:  Entropy  transfer  through  mechanical  interactions  from  phase  [sub¬ 
script]  to  phase  [superscript]  per  unit  mass  [J/(kg-s-°K)[ 

paj :  Partial  mass  density  of  j^1  component  of  ct-phase  [Kg/m3]  so  that  ea paj 
is  the  total  mass  of  j^1  constituent  in  phase  a  divided  by  the  volume  of 
REV 

pa:  Mass  density  of  a-phase  averaged  over  a-phase  [Kg/m3] 

da‘ :  Induced  curl  of  electric  field  of  species  j  due  to  presence  of  other  species 
within  the  same  phase.  See  Faraday’s  law.  [V/m] 

<r/' :  Effect  constituent  j  of  phase  f3  has  on  the  curl  of  the  electric  held  of 
phase  a  [V/m] 


Appendix  B.  Definition  of  Macroscopic  Bulk  Variables 

The  relationships  between  the  macroscopic  constituent  variables  and  their 
bulk  counterparts  follow: 

N 

Pa  =  (R1) 
3= 1 


Caj 


(B.2) 


Swelling  Porous  Media  with  Electroquasistatics 


28 


va 

N 

3= 1 

N 

(B.3) 

■ — -O' 
e0 

=  Y.ca%3’ 

3=  1 

N 

(B.4) 

Ja 

=  Ejaj 

i=i 

N 

(B.5) 

Ja 

=  jai  +QeiV°‘i<a, 

3= 1 

N 

(B.6) 

Da 

=  Y.D°3 

3= 1 

N 

(B.7) 

pa 

=  J2paj 

3  =  1 

N 

(B.8) 

Ea 

Wi. 

II 

(B.9) 

€ 

\Ai  * 

ii 

(B.10) 

^a 

dp 

=  £2? 
i=i 

N 

(B.ll) 

^  a 

3  <=a 

<b 

II 

(B.12) 

H" 

=  V '(H“i  -  P“i  X  t)a’'’a) 
i=i 

N 

(B.13) 

^a 

hp 

=  £S? 

j= i 

N 

(B.14) 

za 

=  J2cajzaj 

3  = 1 

N 

(B.15) 

^a 

ze 

=  J2caj(z7  +ex03zai’a) 

3=1 

N 

(B.16) 

ta 

i 

II 

(B.17) 
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=  Y.c"Jynj- 


(B.18) 


=  Y,°aj9? 


£c“'  (f?  +r/v°>-)  +  Pr„E, 

3= 1  ^ 


(B.19) 

(B.20) 


(B.21) 


=  53  Caj  Maj , 


(B.22) 


53  Cai  ( eai  +  -vai’a  ■  vai’a 


(B.23) 


=  53  Qaj  +  taj  ■  vc 


«j,a  _  p<Xiv<Xj,<*  e<*3  +  -(va3’af  +  vai^pai  .  Et  , 

V  2  / 


=  J2cai(haj  +  (gaj+9?)-va^a), 


n  r  /  i  \ 

53  Q°0  +  T?  •  vaj ,a  +  (^ea-a  +  -(va*’afj 


'^■CX.  „  '^'Ot 

-dpEt  ■  va  -  h0  ■  Et 


53^^, 


(B.24) 

(B.25) 


(B.26) 


(B.27) 


(B.28) 


=  53  Caj  baj , 


(B.29) 


53  +  e“VJ’“ 


(B.30) 
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N 

A"  =  ^  C'“J  A0j , 


(B.31) 

(B.32) 


Appendix  C.  Relationship  between  Macroscopic  and  Microscopic 
Variables 

For  notational  brevity,  define  the  volume  average  over  phase  a  with  angle 
brackets: 


<  ^  >a  (x,t)  =  |(^r-|  J§yf  (r^ha(r,t)  dv(£)  (C.l) 

so  that  if  for  example,  pi  is  constant  throughout  phase  a  then  <  pi  >a  would 
be  that  constant  density. 

Similarly  we  define  the  mass  average  to  be 

i>ia(x,t)  =  <  pJ  \§v  |  J  P’(r,t)‘ip3(r,t)'ya(r,t)  dv(£).  (C.2) 

Thus  for  example,  while  the  volume  average  of  velocity  makes  no  physical 
sense,  the  mass  average  (momentum)  does. 

The  relationships  between  the  macroscopic  constituent  variables  and  their 
microscopic  counterparts  follow.  There  are  no  assumptions  made  about  small 
fluctuations  about  the  average.  What  follows  is  exact. 


pa ’ 

=  <^>a 

(C.3) 

Vai 

- rOt 

= 

(C.4) 

If&j 

- rOi 

ip? 

(C.5) 

Ad, 

e0 

=  - — T  [  P’’(wi,n  ~  v •  na  da 

P^\5Va\J6Aa0  y 

(C.6) 

(C.7) 

=  <P>a 

(C.8) 

J,X’> 

=  <Jj>a  +  <  qJevj  >a  -tg>va* 

(C.9) 

Da-> 

=  <  Dj  >a 

(C.10) 

pctj 

=  <  PJ  >a 

(C.ll) 

Eaj 

=  <  Ej  >“ 

(C.12) 

Et 

=  V  £a  <E>a=  A-  f  Edv 

V  \sv\  Jsv 

(C.13) 

qV 

=  <qi>a 

(C.14) 

ts>  P 
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(C.35) 
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